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Abstract 

The late time evolution of Friedmann-Robertson- Walker (FRW) models 
with a perfect fluid matter source is studied in the conformal frame of / (R) 
gravity. We assume that the corresponding scalar field, nonminimally coupled 
to matter, has an arbitrary non-negative potential function V {(f)). We prove 
that equilibria corresponding to non-negative local minima for V are asymp- 
totically stable. We investigate all cases where one of the matter components 
eventually dominates. The results are valid for a large class of non- negative 
potentials without any particular assumptions about the behavior of the po- 
tential at infinity. In particular for a nondegenerate minimum of the potential 
with zero critical value we show that if 7, the parameter of the equation of 
state is larger than one, then there is a transfer of energy from the fluid to 
the scalar field and the later eventually dominates. 



1 Introduction 

The standard inflationary idea requires that there be a period of slow-roll evolution 
of a scalar field (the inflaton) during which its potential energy drives the universe 

*Talk given in "The Invisible Universe", June 29 - July 3, 2009, Paris. 
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in a quasi-exponential expansion. Besides a cosmological constant, a massive scalar 
field (quintessence) provides the simplest mechanism to obtain accelerated expansion 
of the universe within General Relativity. Therefore, scalar fields play a prominent 
role in the construction of cosmological scenarios aiming to describe the evolution 
of the early and the present universe. Since the nature of the scalar field supposed 
to cause accelerated expansion is unknown, it is important to investigate the gen- 
eral properties shared by all FRW models with a scalar field irrespectively of the 
particular choice of the potential. 

In this paper we study the late time evolution of initially expanding flat and 
negatively curved FRW models with a scalar field having an arbitrary bounded 
from below potential function V ((f)). The scalar field is nonminimally coupled to 
ordinary matter described by a barotropic fluid with equation of state 

p = (7 - l)p, < 7 < 2. 

The corresponding case of initially contracting models was studied in pp. Nonmin- 
imally coupling occurs for example, in string theory [2], in scalar-tensor theories of 
gravity [3], in higher order gravity (HOG) theories [I] and in models of chameleon 
gravity [5]. In particular, for HOG theories derived from Lagrangians of the form 

L = f (R) V=9 + 2L m (vl/) , (1) 

it is well known that under a suitable conformal transformation the field equations 
reduce to the Einstein field equations with a scalar field <fi as an additional matter 
source. 

The conformal equivalence can be formally obtained by conformally transforming 
the Lagrangian ([1]) and the resulting action becomes [6], 




Variation of S with respect to g yields the field equations, 

G^ = T^(^,0)+f^(^,^), (2) 
and variation of S with respect to cf) yields the equation of motion of the scalar field, 

Note that the Bianchi identities imply that 

and therefore there is an energy exchange between the scalar field and ordinary 
matter. 
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2 Flat and negatively curved FRW with an arbi- 
trary non-negative potential 

For homogeneous and isotropic spacetimes the field equations (jSJ) reduce to the 
Friedmann equation, 



and the Raychaudhuri equation, 

while the equation of motion of the scalar field (131), becomes 



+ + y' (0) = l_|2 p . ( 6 ) 

The Bianchi identities yield the conservation equation, 

4-37 • 

p + 3n/pH = — ~^~P^ W 

(see for example P, E]). For simplicity we drop the tilde from all quantities. We 
adopt the metric and curvature conventions of [8]. a (t) is the scale factor, an overdot 
denotes differentiation with respect to time t, H = a/a and units have been chosen 
so that c = 1 = 87rG. Here V (0) is the potential energy of the scalar field and 
V = dV/d(f>. 

In the following we consider initially expanding solutions of (J5])-([7|), i.e. H(0) > 0. 
For flat, k = 0, models the state vector of the system ©-(J7]) is ^0, 0, p, Hj , i.e. we 
have a four- dimensional dynamical system subject to the constraint (j3J). Defining 
y := (p and setting (4 — &y)/y/E =: a, we write the autonomous system as 

= V, 

y = -3Hy-V'(<f>) + ap, 

p = -SjpH - apy, (8) 
1 ..2 7 



subject to the constraint 



3H 2 = p+±y 2 + V(4>). (9) 
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Remark 1 The function W defined by 

W(<f>, y, p,H) = H 2 - 1 - {^-y 2 + V{<j>) + pj , 

satisfies 

W = -2HW, 

and therefore, the hypersurface {W = 0}, representing flat, k = 0, cosmologies is 
invariant under the flow of (jHJ). By standard arguments in ordinary differential 
equations theory it follows that also sgn{W) is invariant under the flow of (jHJ). 
We deduce that solutions with W positive, null, or negative, represent scalar field 
cosmologies with k = —1,0, 1 respectively. Similar arguments applied to the third of 
(jHJ) show that if p > at some initial time to, then p(t) > throughout the solution. 
Denoting by 

e=\y 2 + V{<t>), 

the energy density of the scalar field, we see that 

e + p = -3H(y 2 + 1P ), 

which implies that, for expanding models the total energy e+p of the system decreases. 

The equilibria of (jHJ) are given by (0 = 0*, y = 0, p = 0, H = ±-^/V(0*) /3) 
where V 7 (</>*) = 0. Regarding their stability for expanding cosmologies, we note the 
following facts. Critical points of V with negative critical value are not equilibria and 
they rather allow for recollapse of the model. Moreover, nondegenerate maximum 
points (with non negative critical value) for V are unstable, as can be easily seen 
by linearizing system (jHJ) at the corresponding equilibria and verifying the existence 
of at least one eigenvalue with positive real part. It is interesting to study what 
happens near local minima of the potential with non negative critical value, and a 
stability result can be given when k = 0,-1. 

Proposition 2 Let <fi* a strict local minimum for V(4>), possibly degenerate, with 
nonnegative critical value. Then, p* = ((f)*, y* = 0,p* = 0, if* = \J V ^ ) is an 
asymptotically stable equilibrium point for expanding cosmologies in the open spatial 
topologies k = and k = — 1. 

Sketch of the proof. The proof consists in constructing a compact set Q in M 4 
and showing that it is positively invariant. Applying LaSalle's invariance theorem 
to the functions W and (p + e) in Q, it is shown that every trajectory in Q is such 
that HW — > and H(y 2 + jp) — > as t — > +oo, which means y — > 0, p — > 0, and 
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if 2 — \V{(f)) — > 0. Since H is monotone and admits a limit, V{<p) also admits a limit, 
thus the solution approaches the equilibrium point p*. Details in [9]. ■ 
Similar results were proved in [TU] for separately conserved scalar field and perfect 
fluid. 

3 Energy exchange 

In the following we restrict ourselves to the flat, k — 0, case and study the energy 
transfer from the perfect fluid to the scalar field. We are interested to study the 
late time behavior near the equilibrium point (</> = y = 0, p = 0, H = 
which, by the previous Proposition is asymptotically stable. We suppose that the 
initial data in the basin of attraction of this equilibrium are such that the fluid is 
the dominant matter component, i.e. 

Po > eo, 

and we investigate whether there is a time t\ such that 

e(t)>p(t), Vt>ti. (10) 

This question is relevant to the coincidence problem, that is, why dark energy and 
matter appear to have roughly the same energy density today (see jTTJ and references 
therein) . 

If V (0*) > 0, the transition (fTUj) does always happen. In fact, it easily follows 
that if the critical value of the potential is strictly positive, it behaves as an effective 
cosmological constant and the energy of the scalar field tends to this value whereas 
the energy of the fluid tends to zero. We conclude that the relevant case is when 
V(<^*) = and hereafter we will focus on the case when 0* is a nondegenerate 
minimum. Without loss of generality we suppose that = and therefore, the 
general form of the potential studied can be written in a neighborhood of <j) = as 

^(0) = ^ 2 2 + 0(0 3 ), ^>0. (11) 

Integrating ©, we get 

p{t) = ce- a ^a(t)-^. (12) 

From the Proposition in the previous Section we know a-priori that ((f>, <fi, p, Hj — > 
(0, 0, 0, 0) as t — > oo, hence we can write p(t) ~ ca(t) -37 as t — > oo. 

We now turn our attention to the equation of motion of the scalar field ([6]) with 
the potential ffTTj) . namely 

4> + 3H(p + u 2 (p + 0((p 2 ) = ap. (13) 
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This equation can be solved by the Kryloff-Bogoliuboff (KB) approximation |12j . 
We present an outline of the method for the convenience of readers with no previous 
knowledge of the KB approximation. Consider the differential equation 

4> + rjf (0, 0) + w 2 = 0, 0<r/<l. (14) 

If 77 = 0, the solution can be written as 

(t) = —r cos (— ut + x) and (t) = r sin {—ut + \) > 

U 

where r and \ are arbitrary constants. We are looking for a solution of ( Tl4l) which 
resembles to the form of the simple harmonic oscillator, that is, 

(t) = - r (t) cos (-out + x (t)) and (t) = r (t) sin (-ut + %(*)) . (15) 

CO 

Setting 9 (t) = -ut + \ (t) and substituting (JTHJ) in flU} yields 

dr / 1 \ 

'Vf I — r cos9, r sin9 ) sin^, (16) 



dt \co 

d * V f ( 1 a • a\ a 
—— = / -r cos p, r sin t/ cost/. 

dt r \co J 

For rj small, r(t) and x(t) are slowly varying functions of t. Consequently, in a time 
T = 2tt/uj, r(t) and x(t) have not changed appreciably, while 9 (t) = —ut + x (t) 
will increase approximately by — 2n. Therefore, we replace the rhs of f)16p by their 
average values over a range of 2tt of 9, i.e. the amplitude r(t) is regarded as a 
constant in taking the average (this is the essence of the KB approximation). This 
leads to the equations 

^— = — — [ f ( — rcos9,rsm9 J sin9d9, 
dt 2tt Jo \u J 

^ / f ( _ r cost?, r sin 6U cos 9d9. 



dt 2nr J \u 
If we could apply the KB approximation to ffl3|) with 

77/(0,0) =3H<P + 0^ 2 )-ap, (17) 

then we should obtain for the amplitude of 

/ (3ifwrcost? + C'(r 2 cos 2 ^) -ace- (Q/ " )rcose a- 37 ) sin^t?. (18) 
The integral 

e~ {a/uj)r cos 6 sin 9d9, 



dr _ 1 r271 
~dl~ 2tt 
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vanishes and since if is a slow varying function compared to the rapid oscillations 
of the scalar field, it can be considered as constant during one period, T = 2n/u. 
Therefore, (JT51) becomes 

<t = J-Hr + 0{r% (19) 

since the integration of the quadratic terms vanish. Ignoring third-order terms, eq. 
(TT9|) can be integrated to give r(t) = r (t)a~ 3 / 2 , where r (t) is bounded and bounded 
away from zero uniformly as t — > +00. 

Hence we find that for large t the amplitude of the scalar field varies as a~ 3 / 2 . 
Since the amplitude of (j) has by ffTol) the same time dependence as the amplitude of 
</> and in our approximation 

r-u 12 1 2 1.2 

it follows that 

e = a~ 3 . 

Comparing this result with the time dependence of the fluid density p = a -37 , we 
arrive at the following picture for the evolution of the universe. 

If 7 < 1 the energy density p eventually dominates over the energy density of 
the scalar field e and this universe follows the classical Friedmannian evolution. 
For 7 > 1, since p decreases faster than the energy density e of the scalar field, e 
eventually dominates over p. 

A rigorous proof of these results is given in j9]. The main obstruction to apply 
the KB approximation is due to f lT7|) . indeed, the above argument applies if p/r goes 
to zero. Otherwise, the rhs in ( ITS]) is not infinitesimal and so x{t) could in principle 
be comparable to ut. 

Our results could be interesting in investigations of cosmological scenarios in 
which the energy density of the scalar field mimics the background energy density. 
For viable dark energy models, it is necessary that the energy density of the scalar 
field remains insignificant during most of the history of the universe and emerges 
only at late times to account for the current acceleration of the universe. 
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